Macroscopic mechanical oscillators at the quantum limit using optomechanical cooling 
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We discuss how the optomechanical coupling provided by radiation pressure can be used to cool 
macroscopic collective degrees of freedom, as vibrational modes of movable mirrors. Cooling is 
achieved using a phase-sensitive feedback-loop which effectively overdamps the mirrors motion with- 
out increasing the thermal noise. Feedback results able to bring macroscopic objects down to the 
quantum limit. In particular, it is possible to achieve squeezing and entanglement. 
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(N . I. INTRODUCTION 

> : 

Radiation pressure is a classical efïect which is however at the basis of làser cooling and of many techniques for 
| /' | | the manipulation of quantum states of àtoms. Its classical origin however naturally suggests to extend its application 
to a more macroscopic level. In fact, a number of recent papers ^| have shown how radiation pressure can 

be profitably used to cool macroscopic degrees of freedom as, for example, the vibrational modes of a mirror of an 
optical cavity. The possibility to use optomechanical coupling to cool a cavity mirror in conjunction with a phase- 
t— ( ' sensitive feedback loop was first pointed out in Ref . JÏJ . The optomechanical coupling allows one to detect the mirror 
displacement through a homodyne measurement, and then the output photocurrent can be fed-back in such a way 
■ to realize a real-time reduction of the mirror thermal fluctuations. The scheme prop osed in Ref. § roughly amounts 
to a continuous version of the stochastic cooling technique used in accelerators M. In fact, feedback continuously 
"kicks" the mirror in order to put it in its equilibrium position, and for this reason the feedback scheme of Ref. jïj 
has been called "stochastic cooling feedback" in || f?J . However there are important differences with the scheme used 
çvq , in accelerators and its natural extension in the case of atòmic clouds 

A significant improvement has been then achieved with the first experimental realization Q of the feedback cooling 
scheme. The experiment however employed a slightly different feedback scheme, the so called "cold damping" technique 
which was already known in the study of classical electro-mechanical systems PL Cold damping amounts to applying 
a viscous feedback force to the oscillating mirror. In the experiment of Ref. |2], the displacement of the mirror is 
measured with very high sensitivity, and the obtained information is fed back to the mirror via the radiation pressure 
of anothcr, intcnsity-modulatcd, làser beam, incident on the back of the mirror. 

The experiment of Ref. |^| has been performed at room temperature, where the effects of quantum noise are blurred 
by thermal noise and all the results can be well explained in classical terms (see for example ||). However, developing a 
, fully quantum description of the system in the presence of feedback is of fundamental importance, for two main reasons. 
First of all it allows one to establish the conditions under which the effects of quantum noise in optomechanical systems 
become visible and experimentally detectable. For example, Ref. Q has shown that there is an appreciable difference 
between the classical and quantum description of feedback already at liquid He temperatures. Moreover, a completely 
quantum treatment allows one to establish the ultimate limits of the proposed feedback schemes, as for example, the 
possibility to reach ground state cooling of a mechanical, macroscopic dcgrce of freedom. Quantum limits of feedback 
cooling have been already discussed in [[7, 10 1, where the possibility to reach ground state cooling has been shown for 



both schemes. Here we shall review these results and make a detailed comparison of the cooling capabilities of the two 
different schemes. Morever, the present paper will extend the results of [Q for what concerns the possibility to achieve 
new nonclassical effects in the presence of feedback. In particular we shall see that optomechanical cooling allows 
purely quantum effects, as squeezing and entanglement, to come out also in macroscopic mechanical oscillators. The 
experimental realization of these quantum limits in optomechanical systems is extremely difficult, but the feedback 
methods dcscribcd in this paper may be uscful also for microelectromechanical systems, where the search for quantum 



effects in mechanical systems is also very active 1 1 , [12| [L3| . 

The outlinc of the paper is as follows. In Sec. II we describe the model and derive the appropriate quantum 
Langevin equations. In Sec. III we describe both feedback schemes using the quantum Langevin theory developed 
in ttA ' 



15 . In Section IV we analyze the stationary state of the system, and we determine the cooling capabilities of 
both schemes. Section V is devoted to the various nonclassical effects that can be achieved with the feedback scheme 
of Ref. Q, while Section VI is for concluding remar ks. 
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II. THE MODEL 



We shall consider two examples of optomechanical systems, a Fabry-Perot cavity with an oscillating end mirror 
(see Fig. |] for a schematic description), and a ring cavity with two oscillating mirrors (see Fig. ||). These kind of 
optomechanical systems are used in high sensitivity measurements, as the interferometric detection of gravitational 
waves |^| , and in atòmic force microscopes flf7[ . Their ultimate sensitivity is determined by the quantum fluctuations, 
and thcrcfore each movable mirror has to be described as a single quantum mechanical harmònic oscillator with mass 
m and frequency u) m . The mirror motion is the result of the excitation of many vibrational modes, including infernal 
acoustic modes. The description of a movable mirror as a single oscillator is however a good approximation when 
freqüències are limited to a bandwidth including a single mechanical resonance, by using for example a bandpass filter 
in the detection loop p8| . 

The optomechanical coupling between the mirror and the cavity field is realized by the radiation pressure. The 
electromagnètic field exerts a force on a movable mirror which is proportional to the intensity of the field, which, at 
the same time, is phase-shifted by an amount proportional to the mirror displacement from its equilibrium position. 
In the adiabatic limit in which the mirror frequency is much smaller than the cavity free spectral range, one can 
focus on one cavity mode only, say of frequency u c , because photon scattering into other modes can be neglected 
p9| . Moreover, in this adiabatic regime, the generation of photons due to the Casimir effect, and also retardation 
and Doppler eífects are completely negligible |20j . 

Since we shall focus on the quantum and thcrmal noise of the system, we shall neglect all the technical sources of 
noise, i.e., we shall assume that the driving làser is stabilizcd in intensity and frequency. Including these supplementary 
noise sources is however quite straightforward and a detailed calculation of their effect is shown in Ref. [El| . Moreover 
recent experiments have shown that classical làser noise can be made negligible in the relevant frequency range |2^, ^3) . 

The dynamics of an optomechanical system is also influenced by the dissipative interaction with external degrees 
of frecdom. The cavity mode is damped due to the photon leakage through the mirrors which couple the cavity mode 
with the continuum of the outside electromagnètic modes. For simplicity we assume that the movable mirrors have 
perfect reflectivity and that transmission takes place through a "fixed" mirror only. Each mechanical oscillator, which 
may represent not only the center-of-mass degree of freedom of the mirror, but also a torsional degree of freedom as 
in , or an infernal acoustic mode as in p2| , undergoes Brownian motion caused by the uncontrolled coupling with 
other infernal and external modes at thermal equilibrium. 

Let us first analyze the scheme of Fig. |], which can be described by the following Hamiltonian |2(| 

H = hu) c tfb + huj m (P 2 + Q 2 ) - 2hGb%Q + iKE (&V ÍWoí - be iuJot ) , (I) 

where b is the cavity mode annihilation operator with optical frequency üj c , and E describes the coherent input field 
with frequency ujq ~ uj c driving the cavity. The quantity E is related to the input làser power p by E — -*/ pj c /fru>o, 
being 7 C the photon decay rate. Moreover, Q and P are the dimensionless position and momentum operator of the 
movable mirror M. It is [Q,P] = i/2, and G = (uj c /L)y / h/2muj m represents the optomechanical coupling constant, 
with L the equilibrium cavity length. 

The dynamics of the system can be described by the following set of coupled quantum Langevin equations (QLE) 
(in the interaction picture with respect to UoJob^b) 

Q(t) = Lo m P{t), (2) 
P(t) = -uj m Q{t) + W{t)- lm P{t) + Gb\t)b{t), (3) 
6(t) = - (iüj c -icü + y) b(t) + 2iGQ(t)b(t) + E + ^b in (t) , (4) 

where bi n (t) is the input noise operator [p4[ associated with the vacuum fluctuations of the continuum of modes 
outside the cavity, having the following correlation functions 

(hn(t)b in (t')) = {bl(t)b in (t'))=0, (5) 

(bin(t)bÍ(t'))=5(t-t'). (6) 
Furthermore, W(í) is the quantum Langevin force acting on the mirror, with the following correlation function [ p5| , 

1 7n 



where 



<w(í)w(í')> = — ^{Mt-t')+m(t-t')}, (7) 

Zir uj m i* J 



T r {t) = j duj üj cos(cjé) coth [ — ) , (8) 

\2fcsT/ 
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m) = - f 

Jo 



du) lo sin(wí) , (9) 



with T the bath temperature, 7 m the mechanical decay rate, ks the Boltzmann constant, and w the frequency cutoff 
of the reservoir spectrum. Eqs. (Q), (||), and (^|) show the non-Markovian nature of quantum Brownian motion, which 
becomes particularly evident in the low temperature limit ^7|. The symmetric correlation function becomes 
proportional to a Dirac delta function when the high temperature limit fc^T 3> hzu first, and the infinite frequency 
cutoff limit zu — > oo later, are taken. It is only in this limit that the exact QLE (||)-(U) reduce to the Standard ones 
p4j . The quantum Langevin description given by Eqs. (||)-((§) is morè general than that associated with a master 
equation approach, which needs a high temperature assumption p5[ . 

In Standard applications the driving field is very intense so that the system is characterized by a semiclassical 
steady state with the internal cavity mode in a coherent state \8), and a new equilibrium position for the mirror, 
displaced by G\(3\ 2 /ui m with respect to that with no driving field. The steady state amplitude 8 is given by the 
solution of the classical nonlinear equation 8 = E/(j c /2 + iu> c — ioJo + 2iG 2 /uj m \8\ 2 ). In this case, the dynamics is 
well described by linearizing the QLE around the steady state. Redefining with Q(t) and b(t) the quantum 

fluctuations around the classical steady state, introducing the field phase Y(t) — i ur{t) — b(t)) /2 and field amplitude 
X(t) = (b(t) + b'(t)) /2, and choosing the resulting cavity mode detuning A = uj c — lúq + 2G 2 /ui m 8 2 = (by properly 
tuning the driving field frequency wo), the linearized QLEs can be rewritten as 

Q(t) = uj m P(t), (10) 
P(t) = -uj m Q{t)- lm P{t) + 2G8X{t)+W{t), (11) 

Y{t) - -^Y(t)+2G0Q(t) + ^Y in (t), (12) 
X(t) = -|í(í) + ^I»(í), (13) 

where we have introduced the phase input noise Yi n (t) = i (b\ n {t) — bi n (t)\ and the amplitude input noise -Xj„(í) = 
b\ n {t) + b m (t). 

Similar arguments can be used to describe the dynamics of the ring cavity schcmc with two movable mirrors of 
Fig. 2. Assuming for simplicity that the two movable mirrors are identical (i.e., equal mass, frequency and damping), 
and since the two mirrors are perfectly reflecting, one can easily extend the QLEs for the Fabry-Perot cavity (^)-(Q) 
to the optomechanical system of Fig. 2, and get (in the interaction picture with respect to TuJob^b) 

Qi(t) = u m P!(t), (14) 

A(í) = -uj m Qi(t)+W 1 (t)- lm P 1 (t) + GbHt)b(t), (15) 

Qa(*) - u m Pï{t), (16) 

P 2 (t) = -üj m Q 2 (t)+W 2 (t)- lm P2(t)-GbHt)b{t), (17) 

b(t) = - (iuj c - ilüo + y ) b(t) + 2ÍG (Qi(í) - Q a (t)) b(t) + E + ^f c b m (t) , (18) 

where Qj and Pj (J = 1, 2) are the dimensionless position and momentum operators of the movable mirrors Ml and 
M2, with [Qk,Pj] = (i/2)Skj, and G = (uj c /2\/2L)y / h/2muj m is the new coupling constant, (L now represents the 
equilibrium distance between the movable mirrors, as well as the distance between the fixed mirrors, see Fig. 2). The 
operators Wi and W2 are independent Brownian noise operators acting on the mirrors Ml and M2 respectively, each 
of them with a correlation function specified by Eq. (0). 

It is convenient to adopt the center of mass and relative distance coordinates and conjugate momenta, Q± = 
(Qx ±Q 2 )/\/2, and P± = (P x ±P 2 )/y/2. Then, Eqs.(|l·|)-([ï§) become 

Q+(t) - w m P+(í), (19) 

P+(t) = -uj m Q+(t) + W+(t) - j m P+(t) , (20) 

Q'_(í) = w m P_(t), (21) 

P-(í) = -u m Q-(t) + W-(t)-j m P^(t)+V2Gb\t)b(t), (22) 

b{t) = - (íüj c - iu Q + b(t) + 2V2iGQ-{t)b(t) +E + ^/T c b m (t) , (23) 

where W± = (W 1 ± W 2 )/V2 are again independent Brownian noise operators. It is evident from these equations 
that only the relative motion of the two mrrors is coupled to the radiation field, while the center of mass undergoes 
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a simple Brownian motion. Thus, neglecting the latter, the QLEs for the dynamics of the relative motion and the 
cavity mode are identical to those of the single mirror case, Eqs. (^)-(Q), except for the replacements Q — > Q_, 
P -> P_, and G -> V2G. One can again consider an intensely driven cavity mode, and the fluctuations around 
the semiclassical steady state, as discussed above in the Fabry-Perot case. One has a consequently modified steady 
state coherent amplitude /3, and the linearized equations for the relative motion and the cavity mode (assuming again 
zero cavity mode detuning) become identical to those of the Fabry-Perot case, Eqs. (|ïp|)-(13), with the corresponding 
replacements Q — ► Q_, P — ► P_, /3 — > /3, and G — * v2G. 



III. POSITION MEASUREMENT AND FEEDBACK 



As it is shown by Eq. (|12|), when the driving and the cavity fields are resonant, the dynamics is simpler because 
only the phase quadrature Y(t) is affected by the mirror position fluctuations Q(t), while the amplitude quadrature 
X(t) is not. Thcrcforc the mechanical motion of the mirror can be detected by monitoring the phase quadrature Y(t). 
In the ring cavity case, the phase quadrature Y(t) instead reproduces the relative motion dynamics. However, to fix 
the ideas, we shall focus for simplicity onto the Fabry-Perot cavity case only. The mirror position measurement is 
commonly performed in the large cavity bandwidth limit 7 C 3> G/3, ui m , when the cavity mode dynamics adiabatically 
follows that of the movable mirror and it can be eliminated, that is, from Eq. ( |Ï2] ) , 

Y(t)^^Q(t) + ^B, (24) 



and X(t) ~ -^m(í)/-y/7c from Eq. (|T^). The experimentally detected quantity is the output homodyne photocurrent 

010111 

Y out (t) = 2 Vy /^Y(t) - y/riY? n (t) , (25) 

wherc ?7 is the detection cfficiency and Y^ n (t) is a generalized phase input noise, coinciding with the input noise Yi n (t) 
in the case of perfect detection rj = 1, and taking into account the additional noise due to the incfficicnt detection 
in the general case r\ < 1 JÏ5| . This generalized phase input noise can be written in terms of a generalized input 
noise b v (t) as Y^ n (t) = i [6^(í) — 6^(í)]. The quantum noise b v (t) is correlated with the input noise 6, n (í) and it is 
characterized by the following correlation functions pl| 

(b ri (t)b ri (t')} = {bt 1 (t)b ri (t')}=0, (26) 
(b ri (t)bt l (t'))=5(t-t'), (27) 
(M*)^(0) = (b v (t)b\ n (t')) = Vv5(t - t'). (28) 

The output of the homodyne measurement may be used to devise a phase-sensitive feedback loop to control the 
dynamics of the mirror, as it has been done in Ref. Q, and in Refs. ]2|, |j| using cold damping. Let us now see in 
detail how these two feedback schemes modify the quantum dynamics of the mirror. 

In the scheme of Ref. [Q , the feedback loop induces a continuous position shift controlled by the output homodyne 
photocurrent Y out (t). This effect of feedback manifests itself in an additional term in the QLE for a genèric operator 
0(t) given by 

/7c '« 



O fb (t) = / dt'G mf (t')Y out (t - t') [g mf P(t),0(t)} , (29) 

V Jo 

where G m f(t) is the feedback transfer function, and g m f is a feedback gain factor. The implementation of this 
scheme is nontrivial because it is equivalent to add a feedback interaction linear in the mirror momentum, as it could 
be obtained with a charged mirror in a homogeneous magnètic field. For this reason here we shall refer to it as 
"momentum feedback" (see, however, the recent parametric cooling scheme demonstrated in Ref. Q, showing some 
similarity with the feedback scheme of Ref. Q). 

Feedback is characterized by a delay time which is essentially determined by the electrònics and is always much 
smaller than the typical timescale of the mirror dynamics. It is therefore common to consider the zero delay-time 
limit G rn f(t) ~ < 5(t) , which is quite delicate in general However, for linearized systems, the limit can be taken 

directly in Eq. ( p9| ) M, so to get the following QLE in the presence of feedback 

Q(t) = u m P{t) + 9míloY{t) - i^MY^(t) , (30) 
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P(t) = -w m Q{t)- lm P{t) + 2G(3X{t)+W{t), (31) 



Y(t) = -±Y(t)+2GpQ(t) + ^Y in (t), (32) 

X(t) = -^W + y^W, (33) 

where we have used Eq. (|2^). After the adiabatic elimination of the radiation mode (see Eq. (pi))), and introducing 
the rescaled, dimensionless, input power of the driving làser 



16G 2 /3 2 64G 2 

< = = T jP> (34) 

7 m 7 c huolmlé 



and the rescaled feedback gain <?i = —AGj3g m f /jm, the above equations reduce to 

Q(t) = u m P(t) - lm9l Q{t) - J^- 9l Y in {t) + J^ 9 -±Y? n (t) , (35) 

1 



P(í) = -u; m Q{t)- lm P(t) + -yff^CX m (t)+W(t). (36) 

This treatment explicitly includes the limitations due to the quantum efficiency of the detection, but neglects other 
possible technical imperfections of the feedback loop, as for example the electrònic noise of the feedback loop, whosc 
effects have been discussed in . 

Cold damping techniques have been applied in classical electromechanical systems for many years 0, and only 
recently they have been proposed to improve cooling and sensitivity at the quantum level ||29fl . This technique 
is based on the application of a negative derivative feedback, which increases the damping of the system without 
correspondingly increasing the thermal noise |], |2^] . This technique has been succesfully applied for the first time to 
an optomechanical system composed of a high-fmesse cavity with a movable mirror in the experiments of Refs. ^ |3|, Q. 
In these experiments, the displacement of the mirror is measured with very high sensitivity |2^], and the obtained 
information is fed back to the mirror via the radiation pressure of another, intensity-modulated, làser beam, incident 
on the back of the mirror. Cold damping is obtained by modulating with the time derivative of the homodyne signal, 
in such a way that the radiation pressure force is proportional to the mirror velocity. The results of Refs. g], 
referred to a room temperature experiment, and have been explained using a classical description. The quantum 
description of cold damping has been instead presented in p9[ ] using quantum network theory, and in ^ |7j using a 
quantum Langevin description. In this lattcr treatment, cold damping implies the following additional term in the 
QLE for a genèric operator 0(t), 

i f* 

O fb (t) = — — / dt'G cd (t')Y out (t - t') [g cd Q(t),0{t)} . (37) 
V\Tk Jo 

As in the previous case, one usually assume a Markovian feedback loop with negligible delay. Since one needs a 
derivative feedback, this would ideally imply G cc j(í) = —5'(t), i.e., G c( í(lü) = ílü, Vcj, even though, in practice, it is 
sufhcient to satisfy this condition within the detection bandwidth Au. In this case, the QLEs for the cold damping 
feedback scheme become 

Q(t) = co m P(t), (38) 

P(t) = - UJm Q(t)-j m P(t) + 2GpX(t)-g cd Y(t) + -^ = Yl(t)+W(t), (39) 



Y(t) = -^Y(t) + 2Gf3Q(t) + ^Y m (t), (40) 

X(t) = -^ X (t) + ^X in (t). (41) 

Adiabatically eliminating the cavity mode, and introducing the rescaled, dimensionless feedback gain 172 = 
4G/3w„ l 3 C(i /7„ l 7 c , one has 

Q(t) = Lü m P(t), (42) 



P{t) = -w m Q{t) ~ l m P{t) + - V7mCX in (í) + W(t) 

-^Q(Í) - 9 -^Y in (t) + f^=Yl{t). (43) 
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The presence of an ideal derivative feedback implies the introduction of two new quantum input noises, Yi n (t) and 
Y?L(t), whose correlation functions can be simply obtained by differentiating the corresponding correlation functions 
of Yi n (t) and YP n (t). We have therefore 



(Y m (t)Y m (t')) = (Y in (t')Y in (t)) = (YZ(t)Y&(l>)) = (Y t l{t')Yl(t)) = -5(t t% 

(Yimn(f)) = (Y in (t')Ym) = -vvht-n 

{X in {t)Yl{t')) = -(Yl(t')X in (t))=-iVvS(t-t'), 



(44) 
(45) 
(46) 



which however, as discussed above, have to be considered as approximate expressions vàlid within the detection 
bandwidth only. 

The two sets of QLE for the mirror Heisenberg operators, Eqs. (|35|)-(|36|) and (^2|)-(p3|), show that the two feedback 
schemes are not exactly equivalent. They are however physically analogous, as it can be seen, for example, by looking 
at the differential equation for the displacement operator Q(t). In fact, from Eqs. (p5|) and (pq) one gets 



Q(t) + {l+ gi ) lm Q(t) + (ujf yi 

(Ym+lmYZit) 

for the momentum feedback scheme, while from Eqs. (E2h and (ES) one gets 




gi \ Yin(t) + JmYin{t) 



(47) 



Q{t) + {l + g 2 ) lm Q{t)+Lü 2 m Q{t) 



(48) 



Y m (t) 



W7n 



for the cold damping scheme. The comparison shows that in both schemes the main effect of feedback is the modifica- 
tion of mechanical damping 7 m — > 7 m (l + gi) (i — 1,2). In the momentum feedback scheme one has also a frequency 
renormalization uj^ n —* tüf n + 7 r 2 n .9i, which is however negligible when the mechanical quality factor Q = u) m /'y m is 
large. Moreover, the noise terms are similar, but not identical. In particular, the feedback scheme of Ref. Q is also 
subject to the phase noises Y in (t) and Y^(t), while cold damping is not. However, the comparison shows that also 
momentum feedback provides a cold damping effect of increased damping without an increased temperature. 

IV. COOLING AND STATIONARY STATE 

We now study the stationary state of the movable mirror in the presence of both feedback schemes, which is obtained 
by considering the dynamics in the asymptotic limit t — » oo. We shall see that, in both cases, ground state cooling 
can be achieved. 



A. Cold damping 



Now we characterize the stationary state of the mirror in the presence of cold damping. This stationary state has 
been already studied using classical arguments inlfl ^| , while the discussion of the cooling limits of cold damping in 
the quantum case has been recently presented in jïo|j. The results of (Ï(J have been generalized to the case of nonideal 
quantum cfhciency r\ < 1 in Q , and here we shall review these results and compare in detail the cooling capabilities 
of the two feedback schemes. 

The stationary variances can be expressed in terms of the Fourier transforms of the noise correlation functions, and 
using the correlation functions (||), (ph, ÏÏÏh, (J2q)-(28 ), and (E3)-(|6|), one can write 



(Q 2 } st 



= lr> 



du 2 

7T \Xcd{u)\ 
Z7T 



C | g 2 2 IG^HI 2 

4 



2ui, 



■ coth teJ e H 



4??C 

7,H M 



(49) 
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(P 2 )st 



[°° deu uj 2 

/ o T \Xcd{uj) 



2u> r , 



• coth 



Tou 
2k B T 



Ç g | \G cd {u>)Y 
4 477C 



Q[-OT,07] M 



(50) 



where 



w m - ^ 2 + Í^Tm (1 + .92) 



(51) 



is the frequeney-dependent susceptibility of the mirror in the cold damping feedback scheme, and 0r_ W)CT i(a;) is a 
gate function equal to 1 for \üj\ < m and zero otherwise. 

In general, each steady state variance has three contributions: i) the back action of the radiation pressure, pro- 
portional to the input power £; ii) the feedback-induced noise term proportional to the feedback gain squared, and 
inversely proportional to the input power; iii) the thermal noise term due to the mirror Brownian motion. The 
back-action contribution can be evaluated in a straightforward way for both variances, and it is given by 



(Q 2 ) b st = (i* 



I st 
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92) 



(52) 



The feedback-induced contribution generally depends upon the form of |G C£ ;(w)| 2 , where G c d(uj) is a causal function 
(i.e., it is analytic for Imw < 0) and approximatively equal to ílu within the detection bandwidth Alo. The factor 
IXcd(^)| 2 is highly peaked around the mechanical resonance u> — uj mj with width 7 m (l +52)- Since it is usually 
7m(l +32) < Aw, one can safely approximate |G c( z(ü;)| 2 with its value at the resonance peak, |G c d(u;)| 2 ~ w 2 ^, and 
obtain 



) fb - 

I st — 



) fb - 



9Í 



8r/C(l + 92) 



(53) 



The Brownian motion contribution is generally cumbersome and it has been already exactly evaluated for both 
variances in p6| . However, in optomechanical systems, the condition huj m -C k B T is commonly met, and the classical 
approximation coth(huj /2ksT) ~ 2k B T '/hu> can be made in the expression for (Q 2 ) s t, obtaining 



(Q 



2\BM 



k fí T 



2huj m {l +g 2 )' 



(54) 



In the evaluation of (P 2 )f t AI instead, the classical approximation has to be made with care, because, due to the 
presence of the 1J 1 term, the integral ([so]) has an ultraviolet divergence in the usually considered vj — > 00 limit (see 
also Eq. (pi])). This means that, differently from (Q 2 )f t M , the classical approximation for (P 2 )f t M is vàlid only under 
the stronger condition hm <C k B T [ p6| , and that in the intermediate temperature range hm > fcflT> huj m (which 
may be of interest for optomechanical systems), one has an additional logarithmic correction, so to get 



(P 



2\BM 
st 



k B T 1 

2hoj m 1 + 32 



In 



hi 



hm 
2nkT 



(55) 



However, in the common situation of a high Q mechanical mode, this logarithmic correction can be neglected, the 
dependence on the frequeney cutoff m vanishes, and one finally gets 



<Q 2 ) st = (P 2 )st = 



This fact, together with the fact that 



9Ï 



k fí T 



877C 8 2hu v 



1 



I + .92 



(56) 



(PQ + QP} ( 



hm -(Q(t) 2 ) 
t— >oo dt 



0, 



(57) 



implies that the stationary state in the presence of cold damping is an effective thermal state with a mean excitation 
number (n) = 2(Q 2 ) st ~ 1/2, where (Q 2 ) s t is given by Eq. (|5^). This effective thermal equilibrium state in the presence 
of cold damping has been already pointed out in ^ || , within a classical treatment neglecting both the back-action 
and the feedback-induced terms. The present fully quantum analysis shows that cold damping has two opposite effects 
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on the effective equilibrium temperature of the mechanical mode: on one hand T is reduced by the factor (1 + g^) -1 , 
but, on the other hand, the effective temperature is increased by the additional noisc tcrms. 

Let us now consider the optimal conditions for cooling, and the cooling limits of the cold damping feedback scheme. 
Neglecting the logarithmic correction to (P 2 )f t M , the stationary oscillator energy is given by 



U 8t = 2huj m (Q 2 



flUJr 



4(1 + 52) 



92 , 4fcgT 

— + C H — 



(58) 



This expression coincides with that derived and discussed in pToj j , except for the presence of the homodyne detection 
efhciency i], which was ideally assumed equal to one in E3|. The optimal conditions for cooling can be derived in the 
same way as it has been done in [ [ïo|| . The energy U s t is minimized with respect to £ keeping gi fixed, thereby getting 
Çopt — 32/ \fn- Under these conditions, the stationary oscillator energy becomes 



U st = 



hüJr, 



.92 



2 1 



.92 



1 2k B T 1 
y/rj híú m g 2 



(59) 



showing that, in the ideal limit rj = 1, g 2 — > 00 (and therefore £ ~ 32 ^ 00), cold damping is able to reach the 
quantum limit U s t — hLü m /2, i.e., it is able to cool the mirror to its quantum ground state, as first pointed out in JÏ0| . 

Fig. shows the rescaled steady-state energy 2U s t/huj m versus C plotted for increasing vàlues of 172 (&: 32 = 10, b: 
32 = lCr, c: 32 = 10 5 , d: g 2 — 10 7 ), with kBT/fkü rn — 10 5 and 77 = 0.8. For high gain vàlues, ground state cooling 
can be achieved, even with nonunit homodyne detection efficiency. 



B. Momentum feedback 



Onc can proceed in a similar manner for the evaluation of the steady-state variances in the case of the feedback of 
Ref. H. One has again the back-action, the feedback-induced, and the Brownian motion contributions. Differently 
from cold damping however, the variances now strongly depend on the mechanical quality factor Q and each noise 
contributes differently to {Q 2 ) s t and (P 2 ) st . In the Markovian, zero-delay limit it is G m f(u>) = 1, and the feedback- 
induced noise is delta-correlated as the back-action noise. Using the cor relation functions ©, @, ®, and 



one arrives 



at § 



for the back-action contribution, and 



{Q2) » - wrM^y (60) 

{F)st ~ 8 (l + g 1 )(Q 2 +g l y (biJ 



(Q 2 ) f J = # L 1 + lt gl „ (62) 



&v((l + gl)(Q 2 + 9l) , 
,2 )f b = gï Q 2 
}st 8r,C(l+3i)(Q 2 +3i) 



(63) 



for the feedback-induced contribution. 

Finally, for the Brownian motion contribution, one has a situation analogous to that discussed for the cold damping 
scheme. The generally vàlid expressions for the variances have been already exactly evaluated in pGf . However, in 
the limit hoj m <C k B T, we can make the classical approximation coth(?ïu; /2ksT) ~ 2kBT/huj in the expression for 
(Q 2 )st, obtaining 

(0 2 ) B t M = ksT — (64) 

W ht 2hcü m (l+g2)(Q 2 +giY 1 ' 

Also in this feedback scheme, the Brownian motion contribution to (P 2 ) s t has an ultraviolet divergence in the usually 
considered w — > 00 limit 0. Again, the classical approximation for (P 2 )f t M is vàlid only under the stronger condition 
Tim <C k B T pq |, and, in the intermediate temperature range fiw 3> k B T 3> hu m , one has an additional logarithmic 
correction, i.e., 

(P )bm = 7— : -, 2 r + ln — — . (65) 
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Summing up the three contributions, we arrive at [Q 



(Q 2 ) st 



(P 2 )st 



9Ï 



1 + Q 2 



■9i 



8t?C(1 + 5i)(Q 2 + .9i) 

9Ï Q 2 

87?C(l + 5i)(S 2 + .9i) 

, 7m , ( h w 
H — ■ — m 



C | k B T 
8 '2fíío rn 

8 2hco rn 



Q 2 



(i + gi)(Q 2 + giï 
Q 2 



9Ï 



fji 



(i + gi)(Q 2 + gi) 



2nkT 



(66) 



(67) 



These expressions coincide with the corresponding ones obtained in Q using a Master equation description, except 
for the logarithmic correction for (P 2 ) st , which however, in the case of mirror with a good quality factor Q, is quite 
small, cvcn in the intermediate temperature range hm 3> k B T 3> fko m . 

The feedback scheme of Ref. [[FJ has been introduced just for significantly cooling the cavity mirror. Let us therefore 
study the cooling capabilities of this scheme and compare them with those of the cold damping scheme. The stationary 
oscillator energy U st , neglecting the logarithmic correction of Eq. (|67|)), can be written as 



Us 



hLJ m [(Q 2 ) st + (P 2 ) st ] 



TlLúr, 



g\ (l + 2Q 2 + gl ) 

vC (i + si)(Q 2 + ffi) 



+ C + 



4fc B T \ (.9 2 + 2g 2 + gl ) 
hcj m ) (1 + gi) (Q 2 + gi) 



(68) 



It is evident from Eq. ( pq ) that the effective temperature is decreased only if both Q and g\ are very large. At the 
same time, the additional terms due to the feedback-induced noise and the back-action noise have to remain bounded 
for Q — > oo and g\ — * oo, and this can be obtained by minimizing U s t with respect to C keeping Q and g\ fixed. It 
is possible to check that these additional terms are bounded only for very large Q, that is, if Q/gi — > oo and in this 
case the minimizing rescaled input power is ( op t — gi/ ^/rj- Under these conditions, the steady state oscillator energy 
becomes 



1 2k B T 1 



TlLúr, 



.91 



(69) 



showing that, in the ideal limit i] = 1, gi — > oo, Ç ~ gl — > oo, Q/gi — > oo, also momentum feedback is able to reach 
the quantum limit U st = hu) m /2, i.e., it is able to cool the mirror down to its quantum ground state. The behavior 
of the steady-state energy is shown in Figs. ^ and |^, where U s t (in zero-point energy units hu> m /2) is plotted as a 
function of the rescaled input power Ç In Fig. ^, 2U s t/hui m is plotted for increasing vàlues of g\ (a: g\ = 10, b: 



.91 = 10 J 



91 



10 5 , d: gi = 10 7 ) at fixed Q = 10 7 , and with k B T/hu m = 10 5 and r) = 0.8. The figure shows the 



corresponding increase of the optimal input power minimizing the energy, and that for high gain vàlues, ground state 
cooling can be essentially achieved, even with a nonunit detection efficiency. In Fig. |^, 2U s t/hüj m is instead plotted 
for increasing vàlues of the mechanical quality factor Q (a: Q = 10 3 , b: Q = 10 5 , c: Q — 10 7 ) at fixed g\ = 10 7 . The 
figure clearly shows the importance of Q in momentum feedback and that ground state cooling is achieved only when 
Q is sufficiently large. In the limit Q/g± — > oo, the steady state in the presence of momentum feedback becomes very 
similar to that of cold damping, as it can be checked by comparing the expressions for the steady state variances in 
both cases. Under this limit, the two schemes have comparable cooling capabilities, even though it is evident that 
cold damping is more suitable to reach ground state cooling just because the requirement of a very good mechanical 
quality factor is not needed. The possibility to reach ground state cooling of a macroscopic mirror using the feedback 
scheme of Ref. Q was first pointed out, using an approximate treatment, in |3C| ], where the need of a very large 
mechanical quality factor is under lined. Here we confirm this result using an exact QLE approach. 

The ultimate quantum limit of ground state cooling is achieved in both schemes only if both the input power 
and the feedback gain go to infinity. If instead the input power is kept fixed, the effective temperature does not 
monotonically decrease for increasing feedback gain, but, as it can be easily seen from Eqs. (|58| ) and (|68|), there is 
an optimal feedback gain, giving a minimum steady state energy, generally much greater than the quantum ground 
state energy. The existence of an optimal feedback gain at fixed input power is a consequence of the feedback-induced 
noise term originating from the quantum input noise of the radiation. In a classical treatment neglecting all quantum 
radiation noises, one would have instead erroneously concluded that the oscillator energy can be made arbitrarily 
small, by increasing the feedback gain, and independently of the radiation input power. This is another example of 
the importance of including the radiation quantum noises, showing again that a full quantum treatment is necessary 
to get an exhaustive description of the system dynamics H . 
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The experimental achievement ofground state cooling via feedback is prohibitive with present day technology. For 
example, the experiments of Refs. ^, || have used feedback gains up to 32 — 40 and an input power corresponding 
to Ç ~ 1, and it is certainly dimcult to realize in practice the limit of very large gains and input powers. This is 
not surprising, since this would imply the preparation of a mechanical macroscopic degree of freedom in its quantum 
ground state, which is remarkable. 



V. NONCLASSICAL EFFECTS WITH MOMENTUM FEEDBACK 



The explicit dependence upon the mechanical quality factor Q of the steady state variances makes momentum 
feedback less suitable than cold damping for ground state cooling. However, the presence of an additional, in principle 
tunable, parameter, makes the steady state in the presence of the feedback scheme of Ref. jïj richer, and capable of 
showing interesting and unexpected nonclassical effects at the macroscopic level. 

For example, a peculiar aspect of momentum feedback is its capability of inducing steady-state correlations between 
the position and the momentum of the mirror, i.e., the fact that (QP + PQ) st 7^ 0. This correlation can be evaluated 
solving the QLEs Eqs. © -(§§), and one gets § 

(QP + PQ) st _ fC k B T\ 9l Q g\ Q 



8 2Twj m J (l+5i)(Q 2 +5i) 8íK(l+.9i)(Q 2 +.9i) 

Due to the linearization of the problem (see Eqs. (p~0|)— (|Ï2|) ) , the steady state in the presence of momentum feedback 
is a Gaussian state, which however is never exactly a thermal state, because it is always (Q 2 ) st ^ (P 2 ) st and 
(QP + PQ)st 7^ 0. Its phase space contours are therefore ellipses, rotated by an angle 



1 

— aretan 
2 



(QP + PQ) st 

((Q 2 )st-(P 2 )st) 



with respect to the Q axis. The steady state becomes approximately a thermal state only in the limit of very large Q 
(and Q 2 ^ 31), as it can be seen from Eqs. ([jfjj), ( |67| ) and (j7y). This thermal state approaches the quantum ground 
state of the oscillating mirror when also the feedback gain and the input power become very large. 

A first example of nonclassical behavior is that this Gaussian steady state can become a contractive state, which 
has been shown to be able to break the Standard quantum limit in |3l|j , when (QP + PQ) s t becomes negative, and 
this can be achieved at sufficiently large feedback gain, that is, when gi > rjÇ (£ + 4fcgX '/hu) m ) (see Eq. (f7(ï|)). This is 
shown in Fig. |^. 

A second interesting example is that it is possible to achieve steady state position squeezing, that is, to beat the 
Standard quantum limit (Q 2 ) s t < 1/4. The strategy is similar to that followed for cooling. First of all one has to 
minimize (Q 2 ) s t with respect to the input power Q at fixed g\ and Q, obtaining 



Mn _ giQVi + Q 2 + gi k B r q 2 



W > st 4^(l + ffi)(Q 2 +.9i) 2hu Jm (l + g 1 )(Q 2 +g 1 )- [ ' 

This quantity can become arbitrarily small in the limit of very large feedback gain, and provided that g\ ^> Q 2 . That 
is, differently from cooling, position squeezing is achieved in the limit g± — > 00 (implying Ç — > 00), and there is no 

condition on the mechanical quality factor. Under this limiting conditions, (Q 2 ) s t goes to zero as g 1 1/ ' 2 , and, at the 
same time, (P 2 ) s t diverges as g\^ 2 , so that, in this limit, the steady state in the presence of momentum feedback 
approaches the position eigenstate with Q = 0, that is, the mirror tends to be perfectly localized at its equilibrium 
position. The possibility to beat the Standard quantum limit for the position uncertainty is shown in Fig. ^, where 
(Q 2 )st is plotted versus £ for two different vàlues of the feedback gain, g\ = 10 7 (dotted line), and g\ = 10 9 (full linc), 
with Q = 10 4 , kBT/hiü m = 10 5 , and rj = 0.8. For the higher value of the feedback gain, the Standard quantum limit 
(Q 2 )st = 1/4 (dashed line) is beaten in a range of vàlues of the input power Ç 

This capability to achieve quantum squeezing can be exploited to get a third interesting and somewhat unexpected 
nonclassical effect, i.e., the possibility to have an entangled stationary state of the two movable mirrors in the ring 
cavity scheme of Fig. 2. This entangled state corresponds to a squeezed state of the relative motion of the two mirrors, 
in which their positions are strongly correlated. Such a squeezing could be obtained by applying the feedback scheme 
of Ref. Q to the relative motion of the mirror, which amounts to apply a continuous position shift to both mirror 
(see Fig. 2), in such a way that each mirror is shifted exactly opposite to the other. This means having the following 
additional term in the QLE for a genèric operator 0(t) of the two- mirror system, given by 

Ó fb (t) = 1^ f dt'G- f (t')Y out (t - t') \g mf P-(t),0(t)] , (72) 
V Jo 1 J 
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where G m ^(t) and g mf are the corresponding feedback transfer function and gain factor, respectively. Assuming 



again a zero-delay time limit G~ j(í) = S(t), the treatment of the preceding Section for a single movable mirror can 
be easily applied also to the case of this feedback loop acting simultaneously on the two mirrors. In fact, redefining 
the corresponding rescaled, dimensionless, input power of the driving làser 



c 



32G 2 /3 2 



128G 2 

^07m7c 



(73) 



and the corresponding rescaled feedback gain (73 = — 4v / 2G/3g m y/7 m , one gets an expression for the relative distance 
variance {Q 2 _) s t analogous to that of Eq. (|(36|), that is, 



<Q 2 



à 



1 + Q 2 



93 



8í7C( 1 + 53)(S 2 + g 3 ) 



C , k B T 



8 2hWr, 



:i+.93)(q 2 



93) 



(74) 



Entanglement between the two mirrors in the steady state could be demonstrated using one of the sufficient critera 
for entanglement recently appeared in the literaturc Jp3, T 
satisfy is the so-called product criterion fl32| , [35[ |36| , p7|, whic : 
written as 



5l 36, 37í| . The entanglement criterion easiest to 



1, using the oscillator operators defined above, can be 



£ = 16 {Q 2 _)st (Pl)st < 1 



(75) 



where we have defined the marker of entanglement £. This means that the two movable cavity mirrors become 
entangled when the above product of variances of the relative motion and of the center of mass of the system of two 
mirrors goes below a certain limit. The total momentum is not affected either by the radiation field, or by feedback 



(see Eq. (|72[)), and therefore it is simply given by the Brownian motion contribution (PÍ) 



{PD 



BM: which is 



IBM 



k B T 



in 



2nkT 



(76) 



In order to achieve the best conditions for entanglement one follow the same strategy adopted to get position squeezing 
in the Fabry-Perot cavity case, because one has essentially to squeeze the relative distance. One has first to minimize 
Eq. ( [74] ) with respect to the input power £ at fixed 53 and Q, obtaining the optimal input power 



Çopt — s-, 



93 1 + Q 2 



53 



V 



(77) 



giving exactly expression (71) for (Q 2 -)™ 71 , with g\ — > (73. This quantity can become arbitrarily small in the limit 
of very large feedback gainfimplying a very large input power, see Eq. ([t?])) and fixed Q. Since (P 2 }st does not 
depend upon 33, this implies that the marker of entanglement £ can be reduced below 1. This means the possibility 
to entangle the two mirrors by only means of the feedback action. This is shown in Fig. |[ Notice that in such a case 
an extremely high feedback gain is needed, showing the difficulty in obtaining such cffect. This is in apparent contrast 
with what has been obtained in Ref. [p5| showing that entanglement could be achieved without feedback, and in a 
large temperature range, in the frequency domain, within a narrow bandwidth around the mechanical resonance. The 
key point is that in the present model, the feedback-induced entanglement is achieved for the full stationary state, 
integrated over all freqüències, and is not limited to a small bandwidth in the frequency domain. This is a much 
strongcr cffect, and it is not surprising that it is much more difficult to achieve. 



VI. CONCLUSIONS 



We have studied how quantum feedback schemes can be used in optomechanical system to achieve cooling of vibra- 
tional degrees of freedom of a mirror. We have analysed and compared the momentum feedback scheme introduccd 
in Ref. W , and the cold damping scheme of ||, [| |J, [ï(| . The main effect of feedback is the increase of mechanical 
damping, accompanied by the introduction of a controllable, measurement-induced, noise. The increase of damping 
means reduction of the susceptibility at resonance, and the conseqüent suppression of the resonance peak in the noise 
spectrum. We have then shown the possibility to achieve the ultimate quantum limit of ground state cooling with both 
feedback schemes. Cold damping is more suitable for cooling than momentum feedback because its capabilities are 
not influenced by the mechanical quality factor of the oscillators Q. Instead momentum feedback requires a very large 
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Q, and in this limit it tends to coincide with cold damping. However, for fixed Q and for very large feedback gain (and 
input powers), the steady state in the presence of momentum feedback shows unexpected nonclassical features. In 
fact, it can become a position squeezed, or contractive state. Moreover, applying momentum feedback to the relative 
motion of two movable mirrors of a ring cavity, one can even get an entangled stationary state of the two mirrors, 
again in the limit of very large feedback gains and input powers. 

Both ground state cooling and nonclassical effects are achieved in a parameter region which is extremely difficult to 
achieve with present day technology. This is expected, because it would imply the demonstration of genuinc quantum 
effects for a macroscopic mechanical degree of freedom. 
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FIG. 1: Schematic description of a linear Fabry-Perot cavity with the end oscillating mirror M. The equilibrium cavity length 
is L. A cavity mode is driven by an input làser beam. The output field is subjected to homodyne detection (D). The signal is 
then fed back to the mirror motion (dashed line). 




FIG. 2: Schematic description of a ring cavity with two oscillating mirrors Ml and M2. The equilibrium distance between them 
is L (this is also the distance between the fixed mirrors). A cavity mode is driven by an input làser beam. The output field is 
subjected to homodyne detection (D). The signal is then fed back to the mirrors motion (dashed lines). 
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FIG. 3: Rescaled steady-state energy 2U s t/hum versus the rescaled input power Ç, plotted for difïerent vàlues of g 2 (a: g 2 = 10, 
b: g 2 = 10 3 , c: g 2 = 10 5 , d: g 2 = 10 7 ), with k B T/Twú m = 10 5 and r) = 0.8. The optimal input power correspondingly increases, 
and for high gain vàlues, ground state cooling can be obtained. 
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FIG. 4: Rescaled steady-state energy 2U 3 t/hui m versus the rescaled input power Ç, plotted for different vàlues of gi (a: gi — 10, 
b: gi — 10 3 , c: g\ = 10 5 , d: gi = 10 7 ) at fixed Q = 10 7 , and with kBT/hu m = 10 5 and 77 = 0.8. The optimal input power Ç op t 
correspondingly increases, and for high gain vàlues, ground state cooling can be achieved. 




FIG. 5: Rescaled steady-state energy 2U a t/Tiu! 7n versus £ for increasing vàlues of the mechanical quality factor Q (a: Q— 10 3 , 
b: Q = 10 5 , c: Q = 10 7 ) at fixed gi = 10 7 , and with k B T/huj m = 10 5 and r? = 0.8. 
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FIG. 6: Steady state position-momentum correlation —(QP + PQ) s t versus £ for three vàlues of the feedback gain, gi — 10 5 
(a), gi = 10 6 (b) and gi = 10 7 (c). The other parameters are: Q = 10 4 , k B T/huj m = 10 5 and r) = 0.8. 
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FIG. 7: Steady state position variance {Q 2 ) s t versus Ç for two vàlues of the feedback gain, gi — 10 7 (dotted line), and gi = 10 9 
(full line). The dashed line denotes the Standard quantum limit {Q 2 ) s t = 1/4, while the other parameters are: Q — 10 4 , 
k B T/huj m = 10 5 and rj = 0.8. 




FIG. 8: The marker of entanglement £ is plotted versus the feedback gain gi for three vàlues of the mechanical quality factor 
Q = 10 3 (a), Q = 3 x 10 3 (b) and Q = 10 4 (c). The other parameters are: k B T/hiü m = 10 5 and rj = 0.8. 



